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^ . ABSTRACT 

■ Context. Chemical fractionation in the interstellar medium can create isotopologue abundance ratios that differ by many orders of 
' magnitude from the normal isotopic abundance ratios. By understanding how these differences are reflected in astronomical observa- 
tions under various conditions, we try to obtain a new tool to sensitively probe these conditions. 

' , '■ Aims. Recently, we introduced detailed isotopic chemistry into the KOSMA-r model for photon-dominated regions (PDRs) to give 

■ theoretical predictions for the abundance of the carbon isotopologues as a function of PDR parameters. Combined with radiative 
r ^ ' transfer computations for specific geometries, we estimated the possible intensity ratio of the [C II]/[^'^C II] lines. Here, we compare 

, these predictions with new Herschel observations. 

<~i ■ Methods. We performed Herschel/HIFI observations of the [C II] 158 fj,m line in a number of PDRs. In all sources we observed at 

CLf least two hyperfine components of the [^''C II] transition allowing to determine the [C II]/[^''C II] intensity ratio, after some revision 

of the intrinsic hyperfine ratios. Comparing the intensity ratios with the results from the updated KOSMA-r model, we identify cases 
dominated by chemical fractionation and cases dominated by the optical depth of the main isotopic line. 

Results. An observable enhancement of the [C II]/[^'^C II] intensity ratio due to chemical fractionation depends mostly on geometry 
^ ■ and velocity structure, and less on the gas density and radiation field. It is expected to be at maximum for PDR layers that are some- 

I ^ i ' what shielded from UV radiation, but not hidden by surface layers of optically thick [C II]. In our observations the [C II]/[^'^C II] ratio 

for the integrated line intensity was always dominated by the optical depth of the main isotopic line. However, an enhanced intensity 

■ ratio is found for particular velocity components in a few sources: the red-shifted material in the ultracompact H II region Mon R2, 
^ ■ the wings of the turbulent profile in the Orion Bar, and possibly a blue wing in NGC 7023 . The mapping of the [^'^C II] lines in the 

■ Orion Bar allows to derive a C"'" column density map confirming the temperature stratification of the C"'" layer, in agreement with the 
j^-^ , chemical stratification of the Bar. 

^■f^ . Conclusions. Carbon fractionation can be significant, even in relatively warm PDRs, but a resulting enhanced [C II]/[^'^C II] intensity 

^>f-^ ■ ratio is only observable for special configurations. In most cases, a reduced [C II]/[^'^C II] intensity ratio can be used instead to derive 

• ■ the [C II] optical depth to achieve reliable column density estimates that can be compared with PDR models. The C"'" column densities 

for all sources show that at the position of the [C II] peak emission, a dominant fraction of the gas-phase carbon is in the form of C"*". 

I Key words. ISM: abundances - ISM: structure - ISM: clouds 

T— I ■ 

>■ 

• ^ ■ 1 . Introduction of carbon is important even for relatively warm gas, as long 

■ , . . , , . n , , ^ XT , ^ there is a significant fraction of ionized carbon availablqj 
^ ; The isotopic abundance ratios of elements such as C N and O ^.j^^ ^^^^^-^^ ^ ^^^^^^ formation of i^co and depletion of 

. in the mterstellar medium reflect the star formation history of 13^+ temperatures, shifting the fractional abundances of 

the environment, and observations of molecular isotopic line ra- ^j^^ .^^ ^^^^^^^ ^^^.^^^ ^^j^^. ■ ^.^^.^^ 

tios at (sub)mm wavelengths are a common way to measure this .^mpetes with isotope-selective dissociation governed by a dif- 

effect (e.g. Wilson 1999). ferent shield ing through dust, H2 and other CO isotopologues 

Isotopic fractionation of carbon in the interstellar medium, Jvioo<.r ni IBnno'i 

1 i9y-,/i'^y-, • • 1 • 1 u visser et aij i.zuuv). 

i.e the C/ C ratio in various species, IS driven mainly through . , j- • 1 .1 1. 

, J. . . . ±- ' JO ^ systematic discussion of the observational results on car- 

the fractionation reaction u f j ju .1 ,inno^f j 

bon fractionation was provi d ed by Keene et al. ( 1998) for dense 

^^C+ + CO ^ C+ + "CO + 34.8 K (C 1) molecular clouds and lLiszj (1200% with emphasis on more dif- 

... ■ vc » f u- J- 13/-^ . fuse clouds. Through observations of [^^Cl] and ^•^C^^O in the 

providing a significant energy excess for binding the ^■'C to ^ . „ r— — " ■ 1 ,^^ 1 . 

f, 17° — — — —n 1 1 Ao /iK T-u u • 1 f Orion Bar Keene et alj 019980 found little evidence tor chemi- 

the oxygen (ILanger et all 119841) . Thus chemical fractionation , ^ . ' — : : " . — , ,. . , , , 

cal fractionation, but noticed a contradiction between the chem- 

* Herschel is an ES A space observatory with science instruments pro- 

vided by European-led Principal Investigator consortia and with impor- ' All other fractionation reactions for ^''C have lower excess energies, 
tant participation from NASA. therefore being relevant only in cold gas. 
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ical models predicting a slight relative enhancement of "C, 
but a reduction of ^'^C^^O while the observations showed the 
oppos ite behavior. In contrast to the CO results, Sakai et al. 

found clear evidence for fractionation in C2H under 
the dark cloud co nditions of TMC 1 . In line with previous re- 
sults ( Sonnentruc ker et al ■120071: ISheffer et al.ll2007l:fBurgh et alj 
l2007l)lLisztl (i2007i) summarized the clear indications for fraction- 
ation of CO in translucent and diffuse clouds. They confirmed 
that the fractionation reaction ( IC lb dominates the carbon chem- 
istry for densities n^i^ > 100 cm^'' and provided new reaction 
rates for low temperatures, but also showed the importance of 
the CO destruction in reactions with He+ for lower densities. 

In paper I CRollig & Ossenkopf 201 2), we present results 
from the update of the KOSlVIA-r model for photon-dominated 
regions (PDRs) to the full isotopic chemistry, providing isotopo- 
logue ratios for various species. Comparing the fractionation ra- 
tio of the ^•^C+/^'^C+ abundances with the elemental isotopic 
ratio, we note that 

- the fractionation ratio is always higher than or equal to the 
elemental ratio, i.e. '■'C"^ is always underabundant with re- 
spect to C"^. 

- the fractionation ratio equals the elemental ratio at low Ay, 
but increases significantly towards larger Ay- 

These are direct consequences of reaction ( IC lb . No other im- 
portant reaction or isotope-selective photodestruction acts in the 
opposite direction. At low Ay, temperatures are usually high 
enough so that the excess energy plays no role, but once the 
temperature drops below 50 K the back reaction becomes less 
probable and the fractionation ratio increases rapidly. It stays at 
high values as long as enough '■'C'^ ions are available to feed the 
reaction. Deep into the cloud interiors, these ions are formed 
through the dissociative reaction He"^ + '^'CO. The fractiona- 
tion ratio is then controlled by cosmic ray ionization providing a 
roughly constant ^•^C+/^'^C+ enhancement. 

However, two competing processes change the [C Il]/[^'^C ll] 
ratio in PDRs: chemical fractionation raises the [Cll]/[^'^Cll] 
intensity ratio relative to the elemental abundance ratio, while 
the optical depth of the main isotopic line lowers its intensity 
relative to the [^■^C ll] transitions. 

[^•^C II] was first resolved by early heterodyn e observations 
at the Kuiper Airborne Observatory (KAO, e.g. Borei ko et alJ 
Ll988; Boreiko & Betz 1996) showing no increase of the corre- 
sponding [C ll]/[^'^C II] intensity ratio, but a reduction relative to 
the isotopic elemental ratio corresponding to a significant opti- 
cal depth of the main isotopic line reducing the intensity ratio. 
With the low spatial and velocity resolution of these observa- 
tions, it was, however, impossible to determine wether this av- 
erage reduction represents all density or velocity components in 
the beam or whether individual components are dominated by 
either of the two competing processes. Only from additional ob- 
servations with high velocity resolution preferentially covering 
varying conditions it is therefore possible to separate the impact 
of the two processes in PDRs^ 

The HIFI instrument ( de Graauw et al ] l2010l) on-board the 
Herschel satellite (Pilbratt et al. 2010) allowed for the first time 
to systematically observe [^''Cll] with a good sensitivity and 
high spatial resolution in a number of objects. This allows to 
study the ^^C/^'^C ratio in different environments by comparing 
the observations with the theoretical predictions. We detected 
[^■^C II] in four bright PDRs, allowing to look for the signature 
of chemical carbon fractionation in the warm PDR gas or to de- 
termine the optical depth and column density of C+. 



In Sect. 2 we elaborate on the detailed theoretical predictions 
for i2C+/^^C+ and the resulting intensity ratio [C Il]/["C ll] un- 
der varying physical conditions in the frame of the KOSMA-r 
PDR model. In Sect. 3 we present the HIFI observations of [C ll] 
and [^■^C II]. Section 4 discusses the observational results for the 
Orion Bar, Mon R2, NGC 3603, Carina, and NGC 7023 and 
Sect. 5 concludes with a summary of the matches and discrep- 
ancies. 



2. Model predictions for tlie [C ii]/[^^C 11] ratio 

The development of theoretical models for PDRs is an active 
field of research because of the complex interplay of geomet- 
rical constraints, properties of the UV field, and chemical and 
physical processes in the gas (iRollig et al.l l2007l) . Because of 
this complexity, all models with an elaborated description of 
the PDR chemistry are currently restricted to one-dimensional 
geometries, therefore describing PDRs either as plane-parallel 
slabs or as (ensemble of) spheres. However, real sources like the 
prototypical PDR of the Orion Bar, show more complex geome- 
tries, containing elements of both descriptions, e.g. a large-scale 
plane-parallel configuration superimposed on a clumpy small- 
scale structur e, better described by a n ensemble of spherical 
cores (see e.g. lHogerheiide et anil995h . 

We consider both limiting cases in the frame of the KOSMA- 
T PDR model, a code designed to simulate spherical clumps, 
where we can approximate a plane-parallel PDR by considering 
very massive clumps (M > IOOA/q) that form an almost plane- 
parallel structure in their outer layers. This approach allows us to 
use the new implementation of the fractionation network in the 
KOSMA-r PDR model cod e and to easily co mpare the results 
with other PDR models (see lRolhg et al.ll2007h that lack isotope 
networks so far 

2.1. Model setup 

The KOSMA-r PDR model assumes a spherical cloud geome- 
try exposed to isotropic UV illumin ation, meas ured in units of 
the interstellar radiation field, xo (lDrainelll978l) . Model param- 
eters are the hydrogen density at the surface, rismf, UV field, 
X, and cloud mass, M. The density structure assumes a power 
law increase, nnir) — nsuii{r/Ro)^^'^, from the surface to 1/5 
of the outer radius, Rq, and a constant density further in. For 
the results presented here, we keep the metallicity at solar val- 
ues and assume a standard ele mental isotopic ratio of ^^C/^^C of 
67 (IWakelam & Herbsil2008h . The actual elemental isotopic ra- 
tio is subject to systematic variations with galactocentric radius 
(Wilson 1999) but all the sources discussed in Sect. l3.2l have ap- 
proximately solar galactocentric radii. However, even within the 
solar neighborhood there exist noticeable variations with values 
between 57 and 78 where the value of 67 matc hes the ratio ob- 
served in Orion (Lang er & Penziaslll990l ^1993). We therefore 
assume an error bar of up to 15 % for the elemental abundance 
ratio. 

The PDR model iteratively computes the steady-state chem- 
istry, the local cooling and heating processes, and the radiative 
transfer to predict the density of the different species and the 
intensity of the line and continuum radiation emitted by them. 

The restriction of the code to an isotropic illumination is 
not a fundamental limitation when comparing different geome- 
tries. The penetration of an isotropic radiation field can be 
computed from the penetration of a uni-directional radiation 
field by rescaling the optical depth according to the formalism 
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from lFlannervet alJ (Il980l) and vice versa. iRoUig et alJ (l2007h 
showed a very good match of the different PDR models when 
the rescahng was appUed. When re-interpreting the KOSMA-t 
results for a uni-directional illumination, one thus has to stretch 
the width of the outer layers relative to the given scale. 

2.2. The effect ofcfiemical fractionation on the 
[Cu]/[^^Cu] ratio 

Paper I already provides some predictions for clump-integrated 
intensities for ^•^C-bearing species for the special configura- 
tion of a spherical PDRs. To understand in detail how the 
chemical fractionation structure can translate into an observable 
[C ll]/[^'^C II] intensity ratio, we investigate here how the depth- 
dependent carbon fractionation translates into a varying local 
[^■^Cll] and [Cll] emissivity, using the simplified picture of a 
plane-parallel configuration described by two parameters only: 
the density and the impinging UV field. In the resulting one- 
dimensional PDR, different layers reflect a decreasing UV pene- 
tration. This allows us to disentangle the relative contributions of 
these layers to the observable [C Il]/[^'^C ll] intensity ratio. The 
approach ignores all effects resulting from the uncertain details 
of the cloud geometries and of the optical depth of the line along 
the line of sight to the observer in a particular geometry. 

Fig.[T]shows the temperature and abundance structure of the 
dominating carbon bearing species and the resulting [C ll] and 
[^•^C II] line emissivity as a function of the relative radius - or 
depth from the surface into the cloud, respectively - for an ex- 
ample model with a density at cloud surface of 10'' cm~'^ and a 
radiation field of lOOOxo- The local [C ll] line emissivity can be 
computed in terms of a two-level system 



dv — 



with AE = hiy 
i = 2, 



hc?A_ 
87rfczy2 



X N, 



c+ 



ff„exp(-A.E/fcrex) 
gi + Qu exp(-A£;/fcrcx) 



(1) 



k X 91.2K, the statistical weights — 4: and 



IWiese & Fuhl2007h . providing 



edv « 1011 



K kms 



pc 



X Nc+ 



2exp(~91.2K/Tcx) 
l + 2exp(-91.2K/re> 



(2) 



and we obtain the excitation temperature Tex from the de- 
tailed balance includin g collisional excitat ion through 
( Flower & Launavllig? ?*). atomic hydrogen (iLaunav & Rouefij 
Il977h . and electrons (Wilson & Bell 200l) and radiative de- 
excitation and trapping For [^^Cll], this approach sums over 
all hyperfine components. The emissivity represents the inte- 
grated line intensity per hydrogen column density without op- 
tical depth correction. 

In this example model, most of the [Cll] emission stems 
from the cloud surface where all carbon is in ionized form and at 
high temperatures the C+/'^'^C+ fractionation ratio is equal to the 
elemental isotopic abundance ratio so that the [C Il]/['^'^C ll] in- 
tensity ratio also matches the elemental ratio. Deep inside the 
cloud, where most of carbon is in the form of CO and '^C"^ 
is only produced by cosmic ray ionized He+, the intensity ra- 
tio is enhanced by about a factor of ten relative to the elemen- 
tal isotopic ratio, but the C+ abundance is reduced by six or- 
ders of magnitude so that the enhanced ratio is hidden by the 
surface emission. The region where chemical fractionation may 
be detectable is the transition zone from C+ to CO where the 



^ 1011 Kl^ms"^ = 7.11 X 10" 
the [C II] transition 




r/Ro 




ergs cm sr at the frequency of 



r/Ro 



Fig. 1. Top: Abundance distribution of C^ and CO and the 
corresponding '^'^C isotopologues for a model with Usmf = 
10^ cm~'^, 1000 xo and 1000 Mq. The blue curve shows the cor- 
responding gas temperature structure determining the emissivity. 
Bottom: Resulting integrated optically thin line emissivity jedv 
in K kms^V(cm"^ Pc) for the [Cll] and [^^Cll] lines (black 
curves) as a function of the clump radius. The red line shows the 
intensity ratio provided by the two curves. 



C+/'^'^C+ fractionation ratio amounts to a few hundred and the 
intensity is only a factor ten weaker than at the surface (around 
r/Ro — 0.73 in Fig. [T]i- The actual emissivity ratio can be 
boosted to values above the fractionation ratio by radiative exci- 
tation of C+ in the cool and darker inner regions. As C+ "feels" 
the strong line emission from the clump surface, it gains ex- 
citation temperatures above the kinetic temperature while the 
[^■^Cll] excitation temperature closely follows the kinetic tem- 
perature. Due to the exponential term in Eq. |2]this can provide 
up to a factor ten enhancement of the [C Il]/[^^C ll] intensity ra- 
tio for cold gas at moderate densities deeper in the cloud. 

To better evaluate for which conditions we might observa- 
tionally detect chemical fractionation, we compute the differ- 
ence between the [Cll] emissivity and the ['^■^Cll] emissivity 
scaled by the elemental abundance ratio. This quantity does not 
directly reflect the intensity ratio, as shown in Fig. \T\ but en- 
hancements in this difference indicate conditions under which 
an increased C+/'^'^C+ fractionation ratio eventually turns into 
an observable increased [C n]/['^'^C ll] intensity ratio. Figure |2] 
shows the emissivity difference as a function of the distance 
from the outer radius for varying PDR parameters. In the up- 
per panel, we vary the gas density for a constant radiation 
field of lOOOxo while the lower panel shows the impact of dif- 
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Fig. 2. Differences between [Cll] emissivity and [^^Cll] emis- 
sivity multiplied by the standard elemental isotopic ratio of 67 
showing conditions where a increased C+Z^^C"*" fractionation 
ratio turns into an observable increased [Cll]/[^^Cll] intensity 
ratio. The upper plot uses a radiation field of lOOOxo and varies 
the surface gas density: 10'^ cm^^ (black), 10"* cm^'^ (grey), 
10^ cm^"^ (blue), 10^ cm^'^ (red). The lower plot uses a gas den- 
sity of Usuri — 10^ cm^'^ and varies the radiation field: Ixo 
(black), 10x0 (grey), lOOxo (blue), lO^xo (cyan), lO^xo (green), 
lO^Xo (red), lO^xo (orange). The dashed line at 10 K kms ^ 
(cm^-^pc)^^ gives a rough indication what differences might be 
observable with current technology in a deep integration. 



ferent radiation fields for a constant gas density of Usmf — 
10* cm^"^. Whenever the difference is lower than about 1 K 
kms^^ (cm^'^pc)^^ no [Cll] enhancement is practically ob- 
servable, either because the [^'^C ll] emission is too weak or be- 
cause there is no enhanced C+/^'^C+ fractionation ratio. Three 
factors determine the detectability of an enhanced fractionation 
ratio in the optically thin picture. The radiation field has to be 
strong enough to produce a significant column of ionized car- 
bon, the corresponding layer has to be cool enough so that the 
reaction ( IC II ) is effective, and the layer of enhanced C^/^'^C^ 
needs to be thick enough to provide a significant beam filling in 
a finite telescope beam. Paper I has shown that these conditions 
can be met for a wide variety of conditions, provided that there 
is a sufficient mass of dense gas. With at least 1 Mq at densities 
of lO"* cm^"^ a significant fractionation layer occurs even at high 
UV fields. In the upper panel of Figure |2] we see, however, that 
at densities of 10^ cm"-^ or above, the C+-CO transition region, 
producing an enhanced [Cll]/[^'^Cll] intensity ratio associated 



with strong emission, is very thin. The lower panel shows that 
we find only small [C ll] intensities from clumps exposed to ra- 
diation fields of lOxo or below, making every [C ll] observation 
difficult. At all higher radiation fields, fractionation should be 
detectable. When considering the full parameter range, we find 
some degeneracy in the sense that for higher densities fraction- 
ation may also be observable at low UV fields. The ideal en- 
vironment for an increased C+/^'^C+ fractionation ratio is thus 
the mildly shielded gas in classical PDRs with x ^ lOOxo and 
densities in the order of 10'' — 10^ cm"-^. 

2.3. Radiative transfer effects on tlie [Cu]/[^^Cu] ratio 

The [C II] optical depth can be computed in the same way as the 
emissivity (Eq.|2]l through 

1 - exp(~91.2K/rex) 
cm-2 ^^1 + 2 exp(-91.2K/rox) 

(3) 



emissivity (Eq.O through 

frdv = 7.15 X 10-18^^^ 
J cm^^ 

,okms-i 32.9K 

« 7.15 X IQ-is — X iVc+— ^ . 

cm Jcx 

For a Gaussian velocity distribution this is equal to 



-dv = 



Av X 



(4) 



where f denotes the peak optical depth and Av the FWHJVI of 
the velocity distribution. 

Vice versa, we can determine the optical depth from a mea- 
sured [C ll]/[^'^C II] intensity ratio (IR), at any given frequency 
in the line when knowing the C+/^'^C+ fractionation ratio (FR), 
assuming optically thin [^^C ll] lines, and assuming equal exci- 
tation temperatures of [C ll] and [^^C ll]: 



IRu _ 1 - exp(-T[cii]) 

FR T[CII] 



(5) 



For the line integrated intensity ratio, one has to take the broad- 
ening of the [C II] line due to the optical depth into account. For 
a Gaussian velocity distribution and a line-center optical depth 
f < 10, the increase of the line width can be linearly approxi- 
mated by the simple factor 1 + 0.115f, and the average optical 
depth measured through the integrated line intensity i^ 



(t) w 0.64f 



(6) 



To quantify the impact of the line-of-sight optical depth on 
the [C ll]/[^'^C II] intensity ratio we have to switch from the 
semi-infinite plane-parallel description to a finite geometry, con- 
taining a description for the C+ column density, not only in the 
direction towards the illuminating source, but also in the line 
of sight towards the observer. This is naturally taken into ac- 
count when using the clumpy picture that forms the base of the 
KOSMA-T model, i.e. considering the case of spherical PDRs. 
It introduces the clump mass, i.e. the absolute radius of the 
clump, as a third independent parametefl By integrating over 
the full volume of each clump and averaging the intensity over 
the projected area of the spherical model clumps, we mimic 
the situation of a fractal cloud, consisting of an ensemble of 



^ The general, more complex expression, also cov ering the large op- 
tical depth limit, is given e.g. bv lPhillipset"aLl ( fT979l) . 

■* We do not include a separate interclump medium, but rather rep- 
resent interclump material by the equivalent mass of an ensemble of 
small, transient, low-density clumps. 
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Fig.3. Clump integrated intensity ratio / J([C / /([^^C vs. [C ll] intensity from the KOSMA-r model for different 

parameters. The model parameters n,M, and x, are coded as size, shape, and color of the respective symbols. The red horizontal 
line denotes the assumed elemental abundance ratio ['^C]/['^C] of 67. The insert demonstrates the underlying column density ratio 
for the different clumps. For [^'^C ll], we sum over all hyperfine components. 



many clumps, whe re individual structures cannot be resolved 
(IStutzkietal.lll998h . 

In Fig. [3] we compare the ratio of clump averaged integrated 
line intensities, J I{[Cu])dv/ J I{[^^Cu])dv for a wide range 
of densities, clump masses and impinging radiation fields. The 
insert shows the underlying mean column density ratio for the 
clumps producing the observed intensity ratio (see paper I). 

While the column density ratio falls above the elemental 
isotopic ratio for all model clumps, this does not hold for the 
[Cll]/[^'^Cll] intensity ratio. Models with low densities and all 
models with x > 100 show an intensity ratio below the ele- 
mental ratio, down to a value of 38. In these models, the frac- 
tionation ratio in the column densities is close to the elemental 
ratio so that the reduced intensity ratio is a direct measure of the 
average optical depth of the main isotopic line through Eq. |5] 
The reduction of the line integrated intensity ratio by a factor 
two corresponds to a moderate [Cll] optical depth, (r) « 1.6 
or f ~ 2.5, respectively, for the individual clumps. That means 
that in all models exposed to high UV fields an ionized carbon 
layer of about the same column density is formed, independent 
of the other clump parameters. This is naturally explained by the 
transition from C+ to CO at about Ay ~ 0.2. Different PDR 
parameters only change the physical depth and temperature of 
this transition point. Higher [C ll] optical depths only occur if 
multiple clumps add up their optical depth on the line of sight at 
the same velocity without mutual shielding of the impinging UV 
field, i.e. in special geometries. 

In the clump-integrated picture, only models with low UV 
field X and high densities trace the chemical fractionation, show- 
ing an intensity ratio that is enhanced relative to the elemen- 
tal abundance ratio. A detailed comparison shows that even for 
those clumps some line saturation occurs, so that the intensity 
ratio does not measure the true C+Z^'^C^ fractionation ratio. 



but falls somewhat below it. This result is in contrast to the 
pure emissivity considerations in Sect. 12.21 that indicated that 
medium-density clumps produce an enhanced intensity ratio. 

Optical depth effects are always significant in the consid- 
ered geometry with similar column densities of C"*" towards 
the illuminating sources and towards the observer C+ is abun- 
dant down to a shielding gas and dust column of Ay ~ 1 
(see Fig. If all carbon is ionized in this column, it provides 
Nq+ « 2.5 X 10^^ cm~^, corresponding to a [Cll] optical 
depth of 0.8 for a slightly subthermal excitation temperature of 
Tex = 50 K and the line width of 1.5 1cm s^^ used in the model. 
That means that in isotropic configurations with narrow lines, 
as in the clumpy PDR model, [C ll] always turns slightly opti- 
cally thick when the condition for an efficient fractionation is 
met. Hence, we find a relatively constant [C ll] optical depth of 
0.8-1.6 for all clumps that are big enough to contain the criti- 
cal column for Ay — 1 and that are exposed to a radiation field 
strong enough to ionize the carbon. As a consequence, we expect 
a behavior as in Fig. |3]for isotropic PDRs with narrow velocity 
distributions. Most of the fractionation effect in the C+Z^'^C^ 
abundance ratio that was seen in the local [Cll]Z[^'^Cll] emis- 
sivity ratio in Sect. 12.21 is therfore hidden from the observer by 
an optical depth close to unity in this case. 

A recovery of the true fractionation ratio, i.e. the C+Z^^C+ 
column densities, through the observable [Cll]Z[^^Cll] inten- 
sity ratio requires optically thin [C ll] emission, either through a 
much wider velocity dispersion, reducing the [C ll] optical depth 
in a broader line (see Eq. lU, or an anisotropic geometry that 
provides a higher column of gas for the UV shielding than in 
the line-of-sight towards the observer Only in these configura- 
tions we can expect to detect chemical fractionation through the 
[Cll]Z[^^Cll] ratio matching the predictions from Sect. 12.21 In 
contrast we may also find edge-on configurations with a larger 
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column on the line of sight towards the observer than towards the 
illuminating sources (see the Orion Bar in Sect. 14. Il l where the 
[C II] optical depth may grow beyond 2 providing an intensity 
ratio much below the local C+/^'^C+ fractionation ratio. 

Altogether, we must conclude that in most typical PDRs with 
bright [C II] emission, no significant fractionation is expected to 
appear in the C+ column densities, i.e. the fractionation ratio 
will match the elemental abundance ratio while the intensity ra- 
tio is dominated by the optical depth of the [C ll] line pushing it 
below the elemental ratio. In those cases we can only use the in- 
tensity ratio to compute the column density through Eqs. (|5]l 
and (O. This is an essential parameter to compare with mod- 
els for deriving e.g. the ionization degree or the photoelectric 
heating efficiency. Some carbon fractionation may be detectable 
through an increased [C Il]/[^'^C ll] intensity ratio compared to 
the elemental abundance ratio for sources with high densities 
and moderate UV fields in case of a favorable geometry or rela- 
tively broad velocity distributions. However, even in these cases 
it may be impossible to deduce the exact C+/^'^C+ fractionation 
ratio due to the superimposed optical depth effects. 

3. Observations 

3. 1. [^^Cu] spectroscopy 

C+ has a single fine structure transition, [Cll] ^^3/2 — ^ P1/2, 
at 1900.537 GHz, which for ^^C+ is spHt into three hyper- 
fine components due to the unbalanced additional neutron spin. 
The frequenc i es of these transitions have been determined by 
ICooksv et alJ 1198 6') based on the combination of laser mag- 
netic resonance measurem ents with ab-initio computations of 
ISchaefer & Klemml (Il970l) . Our observations from Sect. 14.1.21 
confirm these frequencies with an accuracy of about 3 MHz, as 
the relatively narrow [C ll] and [^"^C ll] lines in the Orion Bar 
match each other and numerous other molecular transitions to 
within 0.5 1cm s^^ in velocity space when using these rest fre- 
quencies. 

We have computed the relative line strengths expected for the 
three hyperfine transitions in local thermodynamic equilibrium 
(LTE), which - for a magnetic dipole transition - are given by 



Table 1. Spectroscopic parameters of the [^■'C ll] P3/2 P1/2 
transition 



S'hfs — 



(2f' + l)(2f +1) 
2/+ 1 



J F I 
F' J' 1 



(7) 



(Eq. 8 in lGarstangl(ll995h . normalized to unity). Here, J' = 3/2 
and J = 1/2 are the total electronic angular momenta for the 
upper and lower states, / = 1/2 is the nuclear spin of the ^■^C 
nucleus, F' and F are the total angular momenta (including nu- 
clear spin) for the upper and lower states, and the curly brack- 
ets denote a Wigner 6j-symbol. This expression - which rests 
on the standard assumptions that magnetic dipole (Ml) radia- 
tion dominates over electric quadrupole (E2) radiation, and that 
the Hamiltonian commutes with all relevant angular momen- 
tum operators - yields a ratio of 0.625 : 0.25 : 0.125 for the 
F'-F = 2- l:l-0:l-l transition^ All spectroscopic 
parameters are summarized in Table [T] 



^ Our line rati os differ from those (0.444 : 0.356 : 0.20) given 
by ICooksv et alj 17986) and used in previous astronomical studies. 
This discr e pancy may account for the anomaly reported recently by 
iGraf etalJ i2012h . who observed aJ' = 2- l:J' = l- 0ra- 
tio in NGC 2024 th at was considerably larger than that predicted by 
ICooksv etal.Ul986h . It a lso requires a minor rev ision to the elemental 
isotopic ratio inferred by iBoreiko & Betzl ( 119961) from observations of 
the F = 2—1 and 1—0 transitions observed toward M42. That study ap- 
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3.2. HIFI measurements 

We p erformed [^^Cll ] obse rvations with the HIFI instru- 
ment (de Graauw et al. 2010) on-board the Herschel satellite 
(Pilbratt et al. 2010) towards four bright PD Rs in the framework 
of th e key proje cts HEXO S ( Bergin & Herschel Hexos Team! 
'201 ll) and WADI ( lOssenkopf et ab,201 L) . The first region is the 
Orion Bar (map center: 5'^35'"20.81", -5°25'17. 1"), an edge- 
on, almost linear PDR at a distance of 415 pc ( Menten et al.l 
2007), formed by material of the Orion Molecular Cloud 1 
(OMC-1) exposed to the radiation from the Trapezium clus- 
ter providing an FUV flux of 4 — 5 x lO^xo in terms of the 
Ibraine (1978) field. Previous, spectrally unresolved o bservations 
of [C I I] in the Orion Bar have bee n carried out by Stace v et"an 
(Il993h and iHerrmann et all (Il997h using the Kuiper Airborne 
Observatory (KAO). In terms of a pure edge-on configuration, 
the Orion Bar PDR is an improbable candidate for the detec- 
tion of ^^C+ fractionation because of the large column density 
of C+ on the line of sight. However, fractionation effects may 
be observable if the [Cll] emission is domin ated by the clumpy 
structure as seen bv lHogerheiide et al] (1 1 9951) and lLis & Schilkd 
(i2003i) . They found clumps with a density well above 10^ cm~^ 
embedded in a more widely distr ibuted " interclump " gas with a 
density of 10"* - 2 x 10^ cm"^ (ISimon el^ll997b . Moreover, 
the gas in front and behind the main Orion Bar should have a 
more favorable geometry. By mapping [^^C ll], we were able to 
distinguish the three regions. 

The second PDR is formed by the molecular interface around 
the ultracompact H II region Mon R2 at a distance of 830 pc 
(6''07"46.2", -6°23'08.0"). Due to the larger distance, the PDR 
is weaker in [C ll] than the Orion Bar, in spite of the larger FUV 
intensity from the central source of about 3 x 10'''xo (Fuent e et al.l 
2010). The almost spherical PDR geometry around the Hll re- 
gion provides a situation close to the plane-parallel model from 
Sect. l2.2l when observing the source center, with two PDR layers 
at high velocities behind and in front of the H II region along one 
line of sight. The v elocity dispersion of up to 10 kms^^ (e.g. 
lGinardetal.ll20I2h promises a low optical depth at high inte- 
grated line intensities making Mon R2 one of the be st candidates 
to detect an enhanced [C Il]/[^^C ll] intensity ratio (I Pilleri et"al] 
2012a). The situation is, however, complicated by the unknown 
local role o f outflows that have been observed on a somewhat 
larger scale (IXu et al.ll2006l) . 

As a third PDR we observed the clu mp MM2 close to 
NGC 3603 at a distance of about 6 kpc (IStoIte et all l2006l) 
(ll''15'"10.89% -61°16'15.2"), illuminated by a F UV intensity 
of ab out lO'^xo from the OB stars in NGC 3603 (iRoUig et al.l 
l20Ilft . The fourth source is the Carina Nebula at a distance of 



plied a correction factor to account for the F = 1 — 1 transition, which 
was not covered by the available instrumental bandwidth. With our re- 
vised line strengths, the correction factor decreases from 5/4 to 8/7, re- 
sulting in a 9.4% increase in the inferred ele mental ratio, from the range 
of 52 — 61 given by Bore iko & Betzl ( 119961) to a range of 57 — 67 . This 
then resolves a small discrepancy, noted bv lBoreiko & Betj ( 11996) . be- 
tween the ratio obtained from observations of [^'^C II] and the value of 
67 ± 3 inferred from observations of ^^CO. 
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about 2.3 kpc (ISmithll2006l) . illuminated by the OB stai's in the 
massive clusters Trumpler 14 and 16 with a FUV flux of about 
4 X lO'^xo (Okadaetal. 2012). Here, we observed two inter- 
faces, the Carina North PDR (10'^43'"35.14% -59°34'04.3") 
close to Trumpler 14, and the Cai'ina South PDR (10'^45™11.54% 
-59°47'34.3") south of Cai-, including the core IRAS-10430- 
5931. In NGC3603 and Carina, the PDRs ai-e known to occur 
at the surface of clumps, partially visible as pillars in the near 
infrared. Therefore, we expect a behavior matching the clumpy 
model from Sect. l2.3l that predicted no observable enhancement 
of the [C ll]/[^^C II] intensity ratio for UV fields above lO^xo- 

The last PDR is the northern filament in NGC 7023 
(2l'^01'°32.4% 68°1 |25.0" ) at a distance of 430 pc 
(Ivan den Ancker et al.l Il997h illuminated by the pre-main- 
sequence B3Ve star HD200775 providing a FUV flux of HOOxo 
at the p osition of our observation ( Joblin et al. 2010; Pill eri et alJ 
l2012bl) . Here, previous observations (e.g Gerin et al.11998 ') show 
a smooth density increase towards the filament and a stratified 
PDR structure so that we may find the situation modelled 
in Sect. 12.21 The [C ll] mapping results have been presented 
already by JobHn et al. (2010). 

All five PDRs were mapped in [C ll] using the on-the-fly 
(OTP) observing mode, however, only the Orion Bar was bright 
enough to detect [^"^C ll] in the map. Therefore, separate deeper 
dual-beam switch (DBS) measurements were performed towards 
the given coordinates in the other four PDRs. Unfortunately, 
the Orion Bar observations used an LO setting that placed the 
F = 1 — 1 transition out of the covered IF range, so that only 
the F = 2 — 1 and F = 1 — components were observed there. 

The data were calibrated with the standard HIFI pipeline in 
the Herschel Common Software System (HCSS Ott 2010) in ver- 
sion 9.0 (2226). All data presented here are given on the scale of 
the antenna temperature corrected for the forward efficiency, . 
The preference for this temperature scale is motivated by the fact 
that all maps show very extended emission. The scale pro- 
vides correct brightness temperatures if the error beam of the 
telescope is uniformly filled by emission of the same magnitude 
as the main beam. The HIFI error beam at 1900 GHz is mainly 
provided by the side lobes of the illumination pattern appear- 
ing at radii of less than two arcminutes (Jellema et al. in prep). 
Scaling to main-beam temperatures, in contrast, would assume 
a negligible contribution of the error beam pickup to the mea- 
sured intensity, increasing the resulting temperatures by a factor 
1.39 (Roelfsema et al. 2012). The two scales will bracket the true 
brightness temperature. As the actual emission decreases with 
growing distance from the sources, the scale will slightly 
overestimate the error beam contribution and underestimate the 
brightness temperature, but due to the very extended emission it 
should be considerably better than the main beam temperature 
scale. We estimate therefore that the antenna temperature under- 
estimates the brightness temperature by 10-20 %. 

The resulting spectra suffer from baseline ripples due 
to instrumental drifts. They were subtracted using the 
Hif iFitFringe task of the HCSS including the baseline op- 
tion with a minimum period of 200 MHz covering all ripples vis- 
ible in the spectra. This method attributes structures wider than 
200 MHz to instrumental artifacts and all narrower structures to 
sky signal. In this way we only get reliable intensity informa- 
tion for lines narrower than 30 kms^^, a condition that is met 
for all our PDRs. In the DBS observations, the resulting spec- 
tra also suffered from self-chopping because the chopper throw 
of 3' often ends at OFF positions that are still contaminated by 
some emission. This affected the detection of the F = 2 — 1 
hyperfine component in NGC 3603 and Carina falling into the 



range of self-chopping features of the main isotopic line. We 
tried to circumvent this problem by looking into the individual 
single-beam switch spectra, that only use one OFF position left 
or right of the source, not combining them. Although they show 
different degrees of contamination from the two OFF positions, 
some contamination always appeared in both spectra, so that we 
obtained no qualitative improvement with respect to a clear de- 
tection of the hyperfine component. Therefore, we show in the 
following only the dual beam switch spectra providing the better 
signal-to-noise ratio. 

With the high signal to noise ratio of the [C ll] line in the 
Orion Bar map we noticed a small zig-zag structure in reported 
positions and the resulting Orion Bar map, originating from the 
alternating direction in the observation of subsequent OTF lines. 
An ad-hoc correction of the pointing information by 1 .4", mu- 
tually shifting the OTF lines relative to each other, resulted in 
a straight Orion Bar structure. This correction will be imple- 
mented in future versions of the HIFI pipeline, but for this paper 
the correction just gives an estimate for the pointing accuracy. 
The original Orion Bar map also suffered from relatively strong 
emission at the OFF position (5^35™44.92% -5°25'17.1") visi- 
ble as absorption feature in numerous points of the map. A sep- 
arate observation of the OFF position was performed using the 
internal HIFI cold load as primary reference and a secondary 
OFF position more than 12' away from the Bar (5'^35™55.0^, 
— 5°13'18.1") as secondary reference. Fortunately, it turned out 
that the baseline from the subtraction of the cold load spectrum 
was good enough so that the secondary OFF position did not 
have to be used because we detected even from this very remote 
position [C ll] emission as bright as 10 K (see F ig. |5]) in agree- 
ment with the large-scale map of iMookeriea et a l. (2003). This 
has to be considered as a very widely distributed emission floor 
in the Orion region. 

4. Observational results 

Table |2] lists the integrated intensities of the [C ll] line and the 
[^■^C II] components for all sources. The relative intensities of the 
[^■^Cll] hyperfine components are roughly consistent with the 
intensity ratios from Table [T] However, in the observations with 
the best signal-to-noise ratio, the Orion Bar peak and NGC 7023, 
we find a somewhat lower F = 2 — 1 /F =1 — ratio of 2 — 2.4 
instead of 2.5. The uncertainty of the integrated intensities con- 
sist of three contributions: radiometric noise, the definition of the 
appropriate velocity integration range, and the uncertainty in the 
baseline subtraction. When measuring the noise in the spectra, 
we found that its contribution to the total uncertainty is actually 
negligible compared to the two other sources of uncertainty for 
all our observations. 

The definition of the best integration range asks for a com- 
promise because broad wings in the [C ll] profile, in particular in 
Mon R2, favor very broad integration ranges to trace the whole 
[C II] emission while the small separation between the [C ll] line 
and the [^■^C ll] F = 2 — 1 transition and wavy baseline struc- 
tures ask for narrow integration ranges. As the wings are not 
detectable in the [^'^Cll] components and their intensity ratio 
forms the main focus of this study, we used common relatively 
narrow integration ranges here (column 2 of Tab.|2]i. In this way, 
we ignore the wing emission that is detectable in the main iso- 
tope line, but compare the same velocity ranges for the isotopic 
ratios. In the extreme case of Mon R2, the integrated [C ll] inten- 
sity is higher by 15 K kms^^, i.e. 4%, compared to the value 
in Tab. |2] when using a very broad integration range including all 
the wings and in the Orion Bar the value would be increased by 



7 



Ossenkopf et al.: [C II] and ["C II] observations in PDRs 



Table 2. Summary of integrated line intensities from the observations and derived parameters. 



source intej 


'ration range 
[kms^M 


[Cii] 
[K kms"^] 


[^^Cii] _F = 2 - 1 
[K kms'^] 


F = 1 - 
[K kms"^] 


F = 1 - 1 
[K kms"^] 




('^[CII]) 

[10^* 


cm-2] 


Orion Bar, peak 


7-13 


857 ±5 


17.2 ± 1.0 


7.8 ±0.5 


— 


30 ±2 


1.9 ±0.2 


10.1 


Orion Bar, ridge'^ 


7-13 


772 ±5 


12.6 ± 1.0 


5.1 ±0.5 


— 


38 ±3 


1.3 ± 0.2 


7.2 


Orion Bar, front"^ 


7-13 


506 ±7 


4.3 ±0.7 


2.1 ±0.6 


— 


69 ± 12 


0.0 ±0.4 


2.3 


Orion Bar, back"^ 


7-13 


529 ±6 


6.3 ± 0.6 


2.9 ± 0.6 




50 ±6 


0.7 ±0.3 


3.9 


Mon R2, total 


5-25 


362 ±5 




2.9 ±0.7 


0.7 ±0.5 


38 ±10 


1.3 ±0.6 


5.7 


Mon R2, blue 


5-12.5 


173 ±3 




2.4 ±0.7 


0.8 ±0.3 


20 ±5 


3.2 ±0.8 


4.8 


Mon R2, red 


12.5 - 25 


188 ±4 




0.7 ±0.5 


-0.3 ±0.4 


170 ± 120 


0.0 ±0.5 


0.9 


NGC 3603 


10-19 


130 ± 2 




1.4 ±0.2 


0.55 ±0.15 


25 ±5 


2.4 ±0.4 


3.6 


Carina N 


-20 - -5 


143 ±4 




1.3 ±0.3 


0.7 ±0.3 


27 ±7 


2.2 ±0.6 


5.0 


Carina S 


-38- -25 


38 ±2 


2.7±0.3-'' 


-0.1 ±0.3 


0.0 ±0.3 


> 22 


< 2.9 


0.23 


NGC 7023 


-1-7 


91 ± 2 


1.07 ±0.05 


0.53 ±0.07 


0.19 ±0.04 


51 ±6 


0.6 ±0.3 


1.0 



Notes. [C II]/[^^C II] intensity ratio. From the sum over all [^''Cll] hyperfine components assuming the ratios from Table[T]to correct for 
non-detected components. Assuming a uniform Gaussian velocity profile for C"*" and ^'^C"'". Average over the « 1000 square-arcsec 
with /([Cll]) > 700 K kms~^ representing the main Bar emission (see contours in Fig.|4l Average over the « 750 square-arcsec with 
450 K kms-^< /([Cll]) < 650 K kms"^ north-west of the Bar. Average over the ~ 900 square-arcssec with 450 K kms"^< J([CII]) < 
650 K kms~^ south-east of the Bar. The identification of this emission as [^"'C II] F = 2 — 1 is questionable because of the small separation 
from the [C II] line. It may rather represent a different velocity component. '^^ Assuming t = 0.0 to compute the column density. 



IRAC 3.6-5.8-8.0/.tfa, [Cll] contours 




7 K kms^^, i.e. 1 %, in this way. For the other sources no wing 
components can be detected due to self-chopping signatures in 
the spectra (see Fig.fTSli. 

The uncertainty that actually determines the error bars given 
in Tab.|2]stems from the baseline subtraction in the spectra. Only 
to a small degree they are given by radiometric noise, but they 
are dominated by correlated residuals in the baselines that are not 
easily described by statistic measures. We performed numerical 
experiments by varying the number of sinusoidal components 
in the Hif iFitFringe task between 1 and 4 and changing 
the size of the window that is masked, i.e. excluded from the 
baseline fit, from the pure integration range to a range that is 
about three times wider in five steps. We inspected all baseline 
fits manually excluding those where self-chopping signatures or 
line wings affected the fit. The uncertainties in Tab. |2] describe 
the variation in the remaining sample in terms of the total cov- 
ered range, i.e. they do not quantify la errors but the extremes 
of that sample. It is obvious that the relative errors are small for 
the Orion Bar with narrow lines while they increase with the 
line width and the resulting larger baseline uncertainty. For the 
figures shown in the next section we used the subjective "best 
baseline fit" from the sample. 

4. 1. Orion Bar 
4.1.1. [Cii] 

[C II] was mapped along a strip perpendicular to the direction 
of the Bar Figure |4] shows the overlay of the integrated in- 
tensity on the corresponding IRAC 8 /im image. [C ll] shows 
a very smooth structure without any indication for dumpiness 
as seen e.g. in the different isotopologues of CO, HCN, HCO"*", 
CS and H2CO bv 'Hogerheiide et al.( (llQgS*) and 'Lis & Schilkg 
(l200l and to a lower degree also in the C91a carbon radio re- 
combination line bv fWvrowski et alJ (Il997h . This indicates that 
[C II] traces mainly the widely distributed, relatively uniform in- 
terclump medium. Position, structure, and overall extent of the 
[C II] Orion Bar emission are similar to the C9 la radio recombi- 
nation line map. That map shows, however, a larger contrast and 
some clumpy substructure within the emission peak. To some 
degree that may be due to insufficient short-spacing of the in- 
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Fig. 4. Integrated [C ll] contours over-laid on the false-color 
IRAC map of the Orion Bai- (red: 8.0 ^im: 0-3000 MJy/sr, green: 
5.8 Aim: 0-7000 MJy/sr, blue: 3.6 fim: 0-10000 MJy/sr). [C ll] is 
mapped in a strip perpendicular to the Orion Bar. Intensities are 
labelled in units of K kms^^. The colors in the IRAC map are 
provided by the signal in the 3.6, 5.8, and 8 /im channels. The 
yellow circle indicates the HIFI beam size. 



terferometer data, but it probably also reflects real differences in 
the excitation of the C"*". As the C91q: emission depends much 
more sens itively on the gas density and temperature than [Cll] 
(see e.g. lNatta et alJ 19941) . the ratio of the two lines may provide 
a way to better characterize the transition between the widely 
distributed gas and the embedded clumps. The [C ll] emission 
peaks approximately 10" deeper in the cloud compared to the 
PAH emission traced by the IRAC 8 /im map. This coincides ap- 
proximately with the p eak of the H2 1-0 S(l) line measured by 
IWalmslev et al.l (I2OOOI) and correlates very well with the C2H 
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Fig. 5. [C II] profile towards the peak of the Orion Bar emis- 
sion (5'>35™20.6% -5°25'5") averaged over a 43" (FWHM) 
Gaussian beam representing the resolution of the previous KAO 
observations (solid line). The dotted spectrum represents the 
widely distributed emission as measured on our secondary OFF 
position 12.3' away from the Orion Bar in a region without 
molecular emission. The dashed lines show Gaussian fits to the 
two profiles. 



emission observed by Ivan der Wiel et al ] (I2009h . The detailed 
comparison of the stratification in the different lines will be the 
topic of a forthcoming paper (Makai et al. in prep). 

We can compar e our data with previo us K AO observations 
ofJC II] in Orion. iBoreiko et all ((l988l) and iBoreiko & Bet3 
(Il996l) obtained frequency-resolved spectra towards 6^C, but 
not towards the Ori on Bar. The spectra l ly unr esolved maps by 
IStacev et aP (Il993h and lHerrmann et al] (1 19971) showed approx- 
imately the same integrated intensity at the pe ak of the Orion 
Bar an d towards O^C. The value from IBoreiko & Betd 
(Il996h of 3.8 X 10-3 erg cm'^ s"! sr'^ is, however, 16% 
lower than the 4.5 x 10^"^ erg cm^^ s^^ sr~^ derived for the 
integra ted [C ll] intensity at the Orion Bar by iHerrmann et al.l 
(Il997h scaling all intensities to match the previous value from 
IStacev et all (1 19931) towards Orion KL. Figure|5]shows our [C ll] 
HIFI profile towards the peak in the Orion Bar after convolution 
to the 43" resolution of the KAO. The integrated intensity of 
585 K kms~^— 4.1 x 10^^ ergcm^'^ s^^ sr^^ falls between the 
two quoted KAO values. Th e discrepancy with the lower value 
from IBoreiko & Bet3 (Il996l) can be partially explained by self- 
chopping of the KAO observations in the very extended emis- 
sion. Our observation of the secondary OFF position (also shown 
in Fig.|5]) gives an estimate for this extended emission providing 
a 6 % correction of the KA O values. The 9 % discrepancy with 
the IHerrmann et all (Il997l) value could be due to the uncertainty 
in the error beam pickup discussed above or the scaling method 
used in their paper. 

As we find almost perfect Gaussian profiles for the lines in 
the Orion Bar (see Fig.|5]) we performed Gaussian fits to all [C ll] 
spectra. In Fig.|6]we show the resulting ma p of ampl itudes and 
widths. As the observations bv lBoreiko & Bet3 (Il996h indicated 
partially optically thick lines, we expect that the lines become 
broader towards the peak of the emission due to optical depth 
broadening. The map, however, does not show any significant 
change of the line width across the Orion Bar The FWHM of 
the line falls between 4.0 and 4.5 kms^^ only turning broader 
for the more diffuse emission in the north. Consequently, the am- 




Fig. 7. Average profile of the two strongest [^'^Cll] hyper-fine 
components in the Orion Bar compared to the [Cll] profile 
scaled by factors of 0.625/67 and 0.25/67. The profiles were 
averaged over all pixels with a [C ll] integrated intensity above 
700K km/s. 



plitude map gives a perfect match to the integrated line map in 
Fig- El As we confirm the significant optical depth in the follow- 
ing section, this indicates that the [C ll] line must be composed 
of many velocity components which are individually optically 
thick, instead of a continuous microturbulent medium that only 
turns optically thick at the largest column density sum. The cor- 
responding map of line velocities only shows a small gradient 
along the Bar with velocities of about 11 kms^^ in the north- 
east and 10.5 kms^^ in the south-west. 



4.1.2. [i-'^Cll] 

Due to the narrow line width in the Orion Bar, we are able to 
detect the F = 2 - 1 and = 1 - lines of [^^C ll] with- 
out noticeable blending with the [Cll] line. Figure [T] shows the 
line profiles of the two [^'^C ll] components when averaged over 
the region with the brightest emission, i.e. with [C ll] intensities 
of more than 700 K kms^^. For a comparison we have added 
the corresponding [Cll] spectra, scaled by factor of 0.625/67 
and 0.25/67. This scaling corresponds approximately to the line 
intensity of the two [^^C ll] transitions that would be expected 
when we assume the canonical abundance ratio and optically 
thin lines. 
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Fig. 8. Contours of the [C ll] line center optical depth, derived 
from the [C Il]/[^^C ll] ratio, overlaid on the integrated [C ll] in- 
tensities observed towards the Orion Bar 




40 20 -20 

A X ["], center at 83.8357, -5.42141 

Fig. 9. Contours of the [C ll] excitation temperature overlaid on 
the derived column densities towards the Orion Bar. 



When considering the line amplitudes, we find that both hy- 
perfine transitions are approximately 2.5 times brighter than ex- 
pected from optically thin [C ll] emission and the normal iso- 
topic ratio. Assuming the canonical abundance ratio, this corre- 
sponds to a line-center optical depth, f of 2.2. The good signal 
to noise ratio of the data allows us to perform Gaussian fits to the 
hyperfine Unes over the full map providing a full map of intensity 
ratios. Using Eq.|5]and assuming no chemical fractionation, we 
can translate this into a map of optical depths. Figure[8]shows the 
resulting line center optical depth, f, overlaid on the integrated 
[C II] intensities. Here, we have added both hyperfine transitions 
to increase the signal to noise. No reliable values can be derived 
for r ^ 0.8 as Eq.|5]turns very flat for low optical depths. Hence, 
all points with lower optical depth have been masked-out. We 
find a good correlation of the optical depth structure with the in- 
tensity map. The map confirms that the area of bright [C ll] emis- 
sion has optical depths above two and the peak optical depth of 
3.4 occurs close to the intensity peak, but we also see a small 
general shift of the optical depth structure relative to the inten- 
sity structure away from the exciting Trapezium cluster This is 
consistent with the picture of hotter gas at the PDR front and 
somewhat cooler gas, but a higher column density deeper into 
the cloud. 

To better quantify this, we use Eq. [3] to compute the exci- 
tation temperatures, and Eq. |2] together with the optical depth 
correction in Eq.|5]to compute the total C+ column density. The 
resulting distribution is shown in Fig. |90 In the range of sig- 
nificant optical depths, the excitation temperatures fall between 



190 and 265 K. For r < 0.8 we assumed a constant excitation 
temperature of 190 K and directly used the optically thin approx- 
imation of Eq. |2]to compute the C+ column density. 

Our excitation temperatures are in agreem ent with the av- 
erage valu e of 220 K and 215 K derived by [Herrmann et al.l 
(Il997h and lSorochenko & TsivilevI ( l2000h . respectively, for the 
Orion Bar Interclump gas at this temperature with a den- 
sity around 10^_cm_^_also explains the C91a o bserva tions 
by Wvrowski et al.l (Il997h . However, lNattaetal.1 (11994 and 
Walmslev et al.1 (l2000l)~ obtain much higher temperatures from 
the radio recombination lines. This could indicate nonther- 
mal emission, dielectronic recombination, or simply different 
gas compon ents c ontributing to the [C II] and C9 1 a emission 
(iNatta et al.lll994tlSorochenko & Tsivilevll2000l) . 

At the scale of the KAO resolution, our approach provides 
a column density of 7 — 8 x 10^* cm^^ being more than two 
times higher than the value of (Herrmann et al.1 d 19971) . consis- 
ten t with our [C III optica l depth exceeding the value derived 
by iBoreiko & Betd (1 19961) by about the same factor. The dis- 
crepancy can be explained to a large degree by their use of the 
line-averaged optical depth, underestimating the peak depth by 
a factor 0.64 for a Gaussian velocity distribution (see Sect. l2.3l l. 
and their use of different spectroscopic parameters, as discussed 
in Sect. 13.11 Our peak C+ column density of 13 x 10^* cm^^ 
implies a configuration where most of the line-of-sight col- 
umn ^ensity_of_the_Orion Bar of Nn w 13 x 10^^ cm^^ (e.g. 
Hog erheiide et al.lll995h forms a PDR with complete ionization 
of the gas phase carbon - abundance X(C) = ?7.(C)/(n(H) + 
2n(H2)) = 1.2 X lO"'' for dense clouds (IWakelam & Herbsll 
2008). A surprising result is, however, the relatively large C+ 
column density in front and behind the Bar. Molecular line 
obser vations (e.g. iHogerheiide et al.l 119951: Ivan der Werf et al.l 
11996^) deduced column densities that are a factor ten lower than 
in the Bar Consequently, at least half of the [C ll] emission in 
the veil must stem from atomic, more diffuse material. 

Fig.|9]also shows the relative shift between the column den- 
sity and the excitation temperature distributions by 5-8" with the 
higher temperatures towards the illuminating source. The new 
observations thus resolve a temperature gradient in C+ consis- 
tent with the overall stratified structure of the Orion Bar PDR. 
The resulting intensity profile then represents a convolution of 
density and temperature structures. 



4.1 .3. Line averages 

To verify the dependence of the results on the assumption 
of the Gaussian velocity distribution, we repeated the analy- 
sis for selected subregions when considering the line-averaged 
[C ll]/[^^C II] intensity ratio, i.e. (t) instead of f . This is also the 
only possible analysis for the other sources where the signal-to- 
noise ratio is lower and the line profiles are highly structured so 
that the assumption of a Gaussian velocity distribution is clearly 
violated. The use of the line averaged optical depths, (r), corre- 
sponds to the assumption of a macro-turbulent situation, where 
[C II] turns independently optically thick in all velocity channels, 
while the f method assumes a microturbulent Gaussian veloc- 
ity distribution for all tracers, where [C ll] turns optically thick 
only in the line center. A preference for one of the methods is 
only possible based on detailed knowlegde of the source struc- 
ture. The (r) results for all sources are included in the last two 
columns in Table |2] 



When scaling the measured intensities to the main beam efficiency, 
i.e. assuming less extended emission, the optical depth does not change. 



but the excitation temperature and C 
30 % relative to the values given here. 



column densities increase by 
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For the Orion Bar, the (r) results always fall below the col- 
umn densities from the line-center optical depth by 30-35 % in 
spite of the almost Gaussian line profiles. A possible cause is the 
systematic difference in the line width between the [^■^C ll] and 
the [C II] lines going beyond the expected optical-depth broad- 
ening. As seen in the average spectrum shown in Fig.|7] but also 
in the map of [^"^C ll] line widths corresponding to Fig. |6] the 
[^■^Cll] line is always 1.0-1.5 kms^^ narrower than the [Cll] 
line. As a consequence, the intensity ratio is even larger than the 
standard elemental isotopic ratio in the line wings (well visible 
in the blue wing in Fig.|7]i. Together with the constant [C ll] line 
width discussed in the previous section, the enhanced intensity 
ratio in the wings indicates different velocity distributions for 
C+ and ^^C+ so that neither the f method nor the (t) method 
can provide accurate results. The 30-35 % of the emission re- 
maining in the wings of the [C ll] line when fitting it with the 
profile of the [^"^C ll] lines provide an estimate for the accuracy 
of the method. 

The (r) values from the average spectra in Table |2] also con- 
firm the spatial gradient of the column densities discussed for 
Fig. |9l When analyzing the bright gas in front and behind the 
Orion Bar by selecting all pixels with a [C ll] intensity between 
450 and 650 K 1cm ^north-west, i.e. in the direction of the il- 
luminating Trapezium stars, and south-east of the Bar, we find 
significantly higher optical depths behind the Bar while the gas 
in front of the Bar is typically hotter by about 20 K, but thinner 
by a factor 1.7. 

None of the average spectra show an indication for an en- 
hanced C"'"/^^C"'' fractionation ratio. This also applies to the 
spectra from the veil in front and behind the Bar where we 
were not able to detect a significant deviation from the elemental 
abundance ratio due to the relatively larger baseline uncertainties 
at lower intensities of the [^''Cll] lines. Our [^•^Cll] observa- 
tions are consistent with negligible carbon fractionation and an 
enhanced optical depth of the main isotopic line for high column 
densities. This is in agreement with a previous detection of the 
[^^Ci] F = 5/2 - 3/2 transition by Keene et al. (1998) in a 
selected clump where the observed ^^C/C ratio showed no frac- 
tionation effect in contrast to the complementary observations 
for i^C^^O and C^^O. 

An indication for an enhanced fractionation ratio is, however, 
seen in the velocity structure when considering the wings of the 
bright line profiles. The broader [C ll] lines showing enhanced 
[C ll]/[^'^C II] intensity ratio in the wings (e.g. in Fig.|7| indicate 
that an additional gas component, with broader velocity disper- 
sion and enhanced fractionation ratio, thus invisible in [^^C ll], 
contributes to the [C ll] profiles. The constancy of the [C ll] line 
width then could be explained by the mutual compensation of the 
relatively larger contribution from the low-velocity dispersion 
Orion Bar material with the optical depth broadening there. To 
sustain the enhanced intensity ratio, the broad velocity gas com- 
ponent should have moderate gas densities of 10'^ — 10"' cm^'^ 
and face moderate UV fields (x ~ 10^, see Sect.|2]). 

4.2. Mon R2 
4.2.1. [Cii] 

The [C II] OTF map of Mon R2 consists of a single strip across 
the ultracompact H II region from north-east to south-west (see 
iFuente et al.ll2010l ). In Fig. [TO]we show the position-velocity di- 
agram for that cut. This is important for understanding the ori- 
gin of the components of the spectral profile in the comparison 
with the [^^C ll] lines. We find a ring-like structure in position- 




d [■■] from 06 07 46.2 -06 23 08.7 

Fig. 10. Position-velocity diagram of the [C ll] line for the mea- 
sured cut through the Mon R2 PDR. 

velocity space, with two well separated components towards the 
center of the H II region. A better understanding of the position 
velocity structure can be obtained by combining t his informa- 
tion with equivalent maps in several other tracers ( Pilleri et~an 
2012a). In particular the comparison with water lines, that show 
up in absorption in the blue component but in emission in the 
red component, suggests that the diagram traces a PDR around 
the H II region that is accelerated by the radiative pressure from 
the central source, so that the back side is red-shifted and the 
front blue shifted relative to the systemic velocity of the medium. 
However, the observations show the existence of some asymme- 
try between the expanding and receding layer and part of the 
central dip could also stem from self-absorption in the optically 
thick line. 

4.2.2. [i^Cii] 

The dual-beam-switch observations of the [^"^Cll] lines were 
performed at the center position of the OTF strip (position 
in Fig. [TOb . i.e. not at the peak of the [C ll] emission. Future 
follow-up observations at the peak position will provide a better 
accuracy. 

Figure [TT] shows the line profiles of the ['^'^C ll] components 
compared to the corresponding [C ll] spectrum, scaled by the 
factor 0.25/67 and 0.125/67 indicating the expected strength of 
the F = 1 — and F = 1 — 1 hyperfine components in case 
of optically thin emission and a canonical abundance ratio (see 
Sect.[TB. Due to the large width of the [C ll] line, the F = 2 - 1 
component falls in the wing of the [C ll] emission, so that we 
cannot estimate its strength. The striking result for the other two 
transitions is the large difference in the spectral shape compared 
to the main isotopic line. The two hyperfine components are con- 
sistent with each other in terms of the spectral shape, peaking at 
about 9 kms"'^ and in terms of the intensity ratio of about a 
factor two. The peak agrees with a flat-top emission part of the 
[C II] line that would be characteristic for optically thick emis- 
sion and is about three times brighter than the scaled blue [C ll] 
component. This is not reflected in the average optical depth (r) 
derived from the integrated line in Table|2]due to the strong mu- 
tual difference in the line shapes. To get an estimate for the two 
main velocity components, we have split the integration interval 
at 12.5 km and include the corresponding results in Table|2] 
This increases the relative uncertainties, but allows to derive in- 
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Fig. 11. Comparison of the profiles of the [^^C ll] hyperfine lines 
in IVIon R2 with the [C ll] profile in the same spectrum scaled by 
the factors 0.25/67 and 0.125/67. The [i^Cii] F=2-l compo- 
nent is blended with the main isotopic line so that is not shown 
here. 



dependent quantitative results for the components that show a 
very different behavior. For the blue-shifted foreground compo- 
nent, we find the highest [C ll] optical depth in our sample. 

The big surprise is the complete absence of redshifted emis- 
sion from the back of the H II region in [^'^C ll]. The red-shifted 
[C II] line has approximately the same brightness as the blue 
component, but any possible [^"^C ll] emission must be weaker 
than expected from optically thin emission and normal abun- 
dance ratios. For both hyperfine components an emission at the 
O.IK level can be excluded in spite of the baseline uncertainties. 
This indicates the direct observation of fractionation in a PDR 
driving an enhanced ^^C+/^'^C+ ratio. However, when consider- 
ing the extreme uncertainties in Tabled stemming from all pos- 
sible baseline subtractions and the possibility that the redshifted 
[C II] line might be optically thin while the blueshifted line has 
an optical depth > 3, the error bars for the ^^C+/^'^C+ fraction- 
ation ratio still cover the standard elemental isotopic ratio so that 
even this case is not unambiguous. 

4.3. NGC 3603, Carina, and NGC 7023 

Figure [T2I shows the spectra for NGC 3603, Carina North, and 
NGC 7023, obtained in the same way as FiglTT]for Mon R2. As 
discussed in Sect.[321 in NGC 3603 and Carina the F = 2 - 1 
component falls into a [C ll] spectral feature resulting from self- 
chopping in the observation, so that it cannot be resolved. The 
two weaker components are resolved. There is no indication for 
a spectral difference with respect to the [C ll] line and the mu- 
tual line ratio reflects the expected value from Sect. 13.11 In both 
sources, the [^"^C ll] lines are about 2.5 times stronger than ex- 
pected from optically thin [C ll] emission and the normal abun- 
dance ratio. The increase indicates an optical depth in the main 
isotopic line of 2.2 — 2.4. 

The equivalent observation of the Carina South PDR, pro- 
vided no detections of [^''C II]. At the frequency of the F = 2 — 1 
transition we find a weak emission component, but we rather at- 
tribute this to a separate [C ll] velocity component in the com- 
plex velocity profile. As the [C ll] line brightness was more than 
three times weaker than in the other sources, the non-detection 
is consistent with the expected line strength in the optically thin 
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Fig. 12. Comparison of the profiles of the [^^C ll] hyperfine 
lines in NGC 3603 MM2 (upper panel), the Carina North 
PDR (central panel), and NGC 7023 North (lower panel) with 
the [C II] profiles from the same positions scaled by the fac- 
tors corresponding to optically thin emission and the standard 
abundance ratio for the resolved hyperfine components. For 
NGC 3603 MM2 and Carina North, the [^^Cii] F=2-l compo- 
nent falls into a self-chopping feature of the main isotopic line 
so that is not shown here. 



case or cases with a moderate [^C ll] line enhancement as in the 
other sources. 

The narrow [C ll] line in NGC 7023 is well separated from 
the [^"^Cll] F = 2 — 1 component, so that this is the only ob- 
servation in which we resolve all three hyperfine components. 
Moreover the narrow lines provide a very small uncertainty for 
the line intensity from the baseline subtraction so that we get 
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quite accurate relative intensities. The average optical depth en- 
hancement of the [^"^C II] lines is moderate, but we also find a 
significantly different behavior of different velocity components 
resolved in the [Cll] line. The [^■^Cll] lines mainly trace the 
component at 2.8 kms~^ corresponding to t he molecular ma- 
terial in the PDR (see e.g. lGerinet al.lll998l) . Here, we find a 
significant optical depth of the [C ll] line in terms of a decreased 
[C ll]/[^'^C II] intensity ratio. The component at 4 1c m was 
identified as separate filament by iFuente et al.l ( 1996h and asso - 
ciated with ionized material in the cavitv bvljoblin et al.l (l2010h . 
It is also seen in [^''C ll] and shows an intensity ratio correspond- 
ing to the elemental isotopic ratio. Finally there is a high ve- 
locity wing with velocities around 1cm s^^. This component 
may show an increased ^^C+/^'^C+ fractionation ratio through 
[^■^C II] intensities below the scaled [C ll] intensity for all three 
hyperfine components. However, our signal-to-noise ratio is not 
good enough to derive a definite statement on the [C Il]/[^'^C ll] 
intensity ratio in that wing. Due to the almost Gaussian veloc- 
ity components, the estimate of the C"*" column density from 
the integrated line parameters should be a reasonable average 
by slightly underestimating the optical depth from the PDR and 
slightly overestimating the contribution from the ionized gas in 
the cavity. For a total column d ensity iVn 1.3 x 10^^ cm^^ at 
the position of our observation jJoblin et alj|2010h this indicates 
that more than half of the carbon in the gas phase is in the form 
ofC+. 



5. Discussion 

The detection of carbon fractionation under dense cloud con- 
ditions remains difficult. For the two species involved in the 
most important carbon fractionation reaction, C+ and CO, the 
lines of the main isotopologue are often optically thick with 
a considerable uncertainty in the optical depth. Consequently, 
it is hard to determine the absolute column density and the 
CO/^'^CO or C+/^"^C+ ratio. Our observations and simulations 
have confirmed this general knowledge for the C+Z^^C"*" ratio. 
Observations of ^■^C^^O are a possible way out, but challeng- 
ing due to the low line intensities, biased towards high column 
densities, and they cannot trace the impact of isotope-selective 
photodissociation provided by the self-shielding of the main iso- 
topologue. 

Our model predictions for PDRs indicate that the observed 
fractionation ratio in all tracers is dominated by the chemi- 
cal fractionation reaction ( IC lb while isotopic-selective photo- 
dissociation plays a minor role. The chemical fractionation in- 
creases the C+Z^'^C"'" ratio relative to the elemental isotopic ra- 
tio in the gas. In terms of the observable [C Il]/[^'^C ll] ratio, this 
could lead to an increase by a factor of a few under favorable 
conditions. They can be provided by inner layers of moderate 
density, moderate UV PDRs or by dense gas that is exposed to 
low UV fields like clumps in an inhomogeneous PDR, embedded 
in and partially shielded by a thinner interclump medium. The 
main condition for an observable increase of the [C Il]/[^'^C ll] 
ratio is, however, that geometry and velocity distribution prevent 
a large [C ll] optical depth. In most geometries, the conditions 
providing an increased isotopic ratio are paired with high opti- 
cal depths of the [C ll] line that rather reduces the [C Il]/[^^C ll] 
intensity ratio below the canonical value. 

Our observations show that most of the [C ll] emission must 
stem from a widely distributed, low-density component in PDRs, 
being more smoothly distributed than many high-density tracers. 
Only in Mon R2, where the H II region is still expanding into 
the surrounding high-density core, a probably relatively homo- 



geneous high-density PDR is formed. The observational results 
on the [C ll]/[^'^C II] intensity ratio roughly follow the predic- 
tions from the models in Sect|2]without providing enough statis- 
tics for a qualitative confirmation yet. For most of the PDRs, 
the average [C Il]/[^'^C ll] intensity ratio is reduced compared to 
the isotopic ratio due to the optical thickness of the [C ll] line. 
Details strongly depend on the exact geometry. Assuming the 
standard elemental ratio for the C+Z^'^C+ abundance ratio, we 
can measure the optical depth and consequently the C+ column 
density through the intensity ratio. The optical depth of the main 
isotopic line varies considerably across the line profile and be- 
tween the different PDRs, showing peak values above three for 
the Orion Bar and Mon R2. 

Indications for an increased [C Il]/[^^C ll] intensity ratio are 
found when studying details of the line profiles allowing to dis- 
tinguish different velocity components in the same source. In 
three cases, we find components showing C+ fractionation in 
a PDR. Unfortunately, none of them allows to unambiguously 
quantify the degree of fractionation due to the large error bars 
from the noise and baseline uncertainties (see Table |2]i. The 
clearest case is the non-detection of [^"^C ll] in the red-shifted 
component in Mon R2. If we interpret the red component as the 
receding inner backside of the expanding H II region, it repre- 
sents a dense, face-on PDR layer which should be at least par- 
tially optically thick in [C ll]. As the front side shows the same 
integrated [C ll] intensity, but no fractionation and a very high 
[Cll] optical depth instead, the actual configuration must be 
asymmetric. The situation if even more puzzling as our mod- 
els predict negligible fractionation for this high-UV field PDR. 
Further high-resolution observations are required to understand 
why the red-shifted material differs from the front material. The 
second case is provided by the wings of the line profiles in the 
Orion Bar, probably created by diffuse gas with a broader ve- 
locity distribution. It may represent extended, possibly diffuse 
H II gas in the Orion veil. To confirm this origin, deeper integra- 
tions in positions apart from the Orion Bar are needed. A final 
case may be given by the blue wing material seen in NGC 7023. 
However, in all three cases the step from a pure detection of frac- 
tionation to a full quantitative assessment requires observations 
with a somewhat better signal-to-noise ratio than obtained here. 
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